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In a recent paper [Phys. Rev. Lett. 97, 070401 (2006)] the transition rate of magnetic spin-flip 
of a neutral two-level atom trapped in the vicinity of a thick superconducting body was studied. 
In the present paper we will extend these considerations to a situation with an atom at various 
distances from a dielectric film. Rates for the corresponding electric dipole-fiip transition will also 
be considered. The rates for these atomic flip transitions can be reduced or enhanced, and in some 
situations they can even be completely suppressed. For a superconducting film or a thin film of a 
perfect conducting material various analytical expressions are derived that reveals the dependence 
of the physical parameters at hand. 
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I. INTRODUCTION 



Atom 



Harnessing the interactions of electromagnetic field and 
matter is one of the ultimate goals of atom optics. One 
promising approach towards control of matter waves on 
small scales is to trap and manipulate cold neutral atoms 
in microtraps near structures microfabricated on a sur- 
face, known as atom chips [l|. Magnetic traps on such 
atom chips are commonly generated either by microfab- 
ricated current-carrying wire or by poled ferromag- 
netic films 0, Q attached to some dielectric or metallic 
body. However, the proximity of the atoms to the surface 
threatens to decohere the quantum state of the atoms 
through electromagnetic field fluctuations. This effect 
arises because the resistivity of the surface is always ac- 
companied by field fluctuations as a consequence of the 
fluctuation-dissipation theorem. For an atom close to the 
surface of a dielectric body these fluctuating fields can be 
strong enough to drive magnetic and electric dipole tran- 
sitions in the atom, as e.g. shown in recent experimental 
studies 0, 0, @] ■ If the atom is in a magnetic or electric 
trap, these flip transitions may lead to atom loss. Such 
transitions are therefore most often undesirable, and we 
want to reduce them or even suppress them completely. 
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FIG. 1: Schematic picture of the setup considered in our cal- 
culations. An atom inside a microtrap is located in vacuum 
at a distance z away from a dielectric slab with thickness 
H . Vacuum is on both sides of the slab. The dielectric slab 
can e.g. be a normal conducting metal or a superconduct- 
ing metal. Upon making a magnetic spin-flip transition or an 
electric dipole-flip transition, the atom becomes more weakly 
trapped and is eventually lost. 



long time (see e.g. Ref. r 



II. THEORY 



In the present paper we intend to explicitly write down 
the flip rate for both of these types of transitions for the 
dielectric slab as shown in Fig. ^ We will e.g. consider a 
normal conducting slab as well as a superconducting slab, 
as described in Ref.Q- To the best of our knowledge, 
this is not done for a general spin or dipole orientation, 
despite the fact that it in principle has been known for a 
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A. Magnetic Spin Transition 

We begin by considering an atom in an initial state \i) 
and trapped at position = (0, 0, z) in vacuum, near 
a dielectric body. The rate of spontaneous and ther- 
mally stimulated magnetic spin-flip transition into a final 
state |/) has e.g. been derived in Ref. Q, 
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where we have introduced the dimensionless components Refs.[Tl|), with the result 



Sj = (f\Sj/h\i) of the electron spin operators Sj corre- 
sponding to the transition \i) — > \f), with j — x,y,z. 
Here \ib is the Bohr magneton, gs ~ 2 is the electron 
spin g factor, and G(ta, ya, is the dyadic Green ten- 
sor of Maxwell's theory. Eq. follows from a consis- 
tent quantum-mechanical treatment of electro mag netic 
radiation in the presence of absorbing bodies |lOt Hl| . 
Thermal excitations of the electromagnetic field modes 
are accounted for by the factor (n th + 1), where n th — 
l/(e hu '^ kBT ~l) is the Planck distribution giving the mean 
number of thermal photons per mode at frequency u> of 
the spin-flip transition. Here T is the temperature of 
the dielectric body, which is assumed to be in thermal 
equilibrium with its surroundings, and ks is Boltzmann's 
constant. The dyadic Green tensor is the unique solution 
to the Helmholtz equation 



V x V x G(r, r', uj) - fc 2 e(r, w)G(r, r', u) = 5(r - r')l . 



(2) 



with appropriate boundary conditions. Here k — lo/c is 
the wavenumber in vacuum, c is the speed of light and 
1 the unit dyad. The tensor G(r, y',oj) contains all rele- 
vant information about the geometry of the material and, 
through the electric permittivity e(r, u>), about its dielec- 
tric properties. Due to causality, any complex dielectric 
function must in general obey the Kramers-Kronig rela- 
tions. Since we only consider non-zero frequencies in a 
suitable finite range, such dispersion relations will be of 
no concern in the present paper. 

The current density in a superconducting media is 
commonly described by the Mattis-Bardeen theory [l2| . 
Following Ref. 0, assuming non-zero frequencies < 
u> <C oj g = 2A(0)/h, where ui is the angular frequency 
and A(0) is the energy gap of the superconductor at 
zero temperature, the current density is well described 
by means of a two-fluid model [131 Il4| . The dielectric 
function is in this case given by |7j 
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Here we have introduced the dimensionless spin fac- 
tor S 2 = S 2 + Sy + S 2 . The unbounded free-space 
lifetime corresponding to this magnetic spin-flip rate is 
ri = l/T°. 

In the following we apply our model to the geometry 
shown in Fig. ^ where an atom is located in vacuum at 
a distance z away from a dielectric slab with thickness 
H . This slab is described by dielectric function as given 
by Eq. ©. The total magnetic transition rate 



r B = (r^ + rjr)(nth + i), 



(6) 



can then be decomposed into a free part and a part purely 
due to the presence of the slab. The latter contribution 
for an arbitrary spin orientation is given by 



F| ab = 2f° + 5 y 2 )/|| + S?^ , (7) 



with the atom-spin orientation dependent integrals 
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The scattering coefficients Cjy(g) and Cm{q) are given 
by 



j _ e i2ri(u>)kH 

C N (q)=r p (q) i _ l2 ^ )kH , 



(10) 



where Xl(T) = e^J m/no n s (T) is the London penetra- and 
tion length and where 5{T) = a n (T) is the skin 

depth associated with the normal conducting electrons. 
As usual, no is the permeability of vacuum and e is the 
elementary charge. The total electron density n (T) is 
constant and given by uq — n s (T) + n n (T) , where n s (T) 
and n n (T) are the electron densities in the superconduct- 
ing and normal state, respectively, at a given temperature 
T. The optical conductivity corresponding to Eq. is 
a(T) = 2/ujfi S 2 (T)+i/ujnoXl(T). Above the transition 
temperature T c , the dielectric function in Eq. © reduces 
to the well known Drudc form. We also stress that the 
theory in this paper is particular to non-magnetic media. 

The rate of a magnetic spin-flip transition for an 
atom in the unbounded free-space is well known (see e.g. 



C M (q) = r s (q) 
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with the electromagnetic field polarization dependent 
Fresnel coefficients 

r*( 9 ) = ^44 , r p ( q ) = f\^-^l (12) 
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Here we have defined rj(uj) = y/ e(uj) — q 2 and 770 = 
y/1 — q 2 . For a thick slab with H = 00, the above equa- 
tions are reduced to the results in Ref. 1151. 



B. Electric Dipole Transition 



III. TWO LIMITING CASES 



The previous section concerns magnetic field fluctua- 
tions. In this section we will consider electric field fluctu- 
ations. For an electrical dipole transition, the rate spon- 
taneous and thermally stimulated decay is given by (see 
e.g. Refs.|I3): 
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where dj 



(f\dj\i), with j = x,y,z, is the matrix ele- 
ment of the atomic dipole operator dj in the direction j, 
corresponding to the transition \i) — ► |/). 

Let us apply this model to the geometry shown in 
Fig. ^ where an atom is located in vacuum at a distance 
z away from a dielectric slab with thickness H . The total 
electric transition rate 



r B = (r°+r s i ab )(n th + i) 



(14) 



can then be decomposed into a free part and a part purely 
due to the presence of the slab. The latter contribution 
for an arbitrary dipole orientation is given by 
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where we have introduced the dipole factor d 2 = d 2 + 
dy + d 2 . The dipole orientation dependent integrals are 
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i.e. the dipole-flip rate in unbounded vacuum for the elec- 
tric dipole transition \i) —> \f), with the corresponding 
free-space lifetime tJ? = 1/rJ,. We mention that Eq. i|18|) 
is consistent with the definition in Refs. [Tl| . 



A. Magnetic Spin Transition 

1. The Limit \ L (T) < S(T),H,\ 

Let us now consider a special case of the dielectric 
function in Eq. J3J. The superconducting term dom- 
inates over the normal conducting term provided that 
Al(T) < 6{T). If, in addition, X L (T) < A, which holds 
true in practically all cases of interest, then we can ne- 
glect the unit term in Eq. ©. Here A = 27r/fc is the 
wavelength associated to the magnetic spin-flip transi- 
tion. The dominant factor in the dielectric function is 
in this case real, and the main contribution to the in- 
tegrals in Eqs. JSJ and © occurs for values of q such 
that < q < 1. For a slab with a thickness such that 
^l(T) >C H, which also holds true in practically all cases 
of interest, the exponential functions in the scattering 
coefficients Eqs. (llOf) and (lllll can be neglected. The 
scattering coefficients are then reduced to Ci\r(g) « r p (q) 
and Cm (q) ~ r s (q) for all relevant values of q. Further- 
more, with the above mentioned assumptions, the Fresnel 
coefficients are reduced to r p (q) w 1 and r s (q) w —1. The 
integrals in Eqs. iJSJ and jSjl can then be solved analyti- 
cally. The total magnetic spin-flip rate for an atom above 
a slab is then 
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where we have defined 
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(21) 
(22) 



Note that Eq. (|20|l is not valid for an arbitrary small 
thickness H of the slab. In the limit Al(T) — > 0, the 
magnetic spin-flip rate in Eq. I|2(J|) is exact. This result 
is consistent with Ref. [l5|. 

Let us now consider the near-field case kz = =?- z <C 1, 
which holds true in practically all cases of interest. The 
magnetic spin-flip rate is then reduced to 



1 B 



{2kzf 
10 



r£(n th + i 



(23) 



provided the atomic spin is oriented perpendicular to the 
slab, i.e. provided that (f\S x \i) — (f\S y \i) = 0. This re- 
sult implies that, for an atom at the surface of the slab, 
i.e. z — 0, there is no magnetic spin-flip at all (see lower 
graph in Fig. [2| . The particular atomic spin orientation 
under consideration is the only one that can give zero 



spin-flip rate despite the presence magnetic field fluctu- 
ations. Furthermore, when the atomic spin is oriented 
parallel to the slab, the magnetic spin-flip rate is 
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2T° (n th + l), 



(24) 



for the near-field case kz <c 1. This result shows that, 
in the near-field regime, the magnetic dipole-flip rate is 
twice the rate as compared to an atom in unbounded 
free-space, as e.g. pointed out in Ref. [3. In current 
atomic chip design (see e.g. Ref.Q) the typical atomic 
frequency is lu/2tt = 560 kHz and a typical atom-surface 
distance is z — 50 /im. Hence, we have kz ~ 10~ 7 , i.e. 
far within the near-field condition fed. 

In passing we also give the small H expansion. For 
sufficiently small thickness of the slab, i.e. H <C 
6 2 (T)/\ L (T), z, Xl(T), the magnetic spin-flip rate is 



t_b ~ r B (n t h + 1, 
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This expression for the magnetic spin-flip rate as it is 
only valid for thicknesses of the slab smaller than the 
London penetration length. 



2. The Limit S(T)<\l(T),H,\,z 

Let us now consider the case when the dielectric func- 
tion as given by Eq. (j^J is dominated by the normal con- 
ducting term, i.e. when 8(T) < X L (T) and 8(T) <C A. 
The dominant factor in the dielectric function is in this 
case imaginary, corresponding to the well known Drude 
form, and the main contribution to the integrals in Eqs. 
(JSJ) and 10 is for values of q such that q < 1/kz. If, in 
addition, 8(T) <C z then the exponential functions in the 
scattering coefficients in Eqs. and are negligi- 
ble for all values q < l/kz provided that 8(T) < H. The 
scattering coefficients are then reduced to Cpf(q) ~ 1 and 
Cm (q) ~ — 1, i-e. the same result as in the last subsec- 
tion. Hence, the conditions 8(T) <C X L (T), H, X, z and 
Xl(T) <C 6(T), H, X give the same results. In particular, 
for the perfect normal conducting limit, i.e. 8(T) — > 0, 
the magnetic spin-flip rate as given by Eq. (|20|l is ex- 
act. Note that the perfect normal conducting limit is 
only valid for the case 8(T) Xl{T),H, X, z, which e.g. 
means that the slab can not be arbitrarily thin. It also 
means that, in contrast to the case as described in last 
subsection, the atom-surface distance z can not be chosen 
arbitrary small. 

We close this section by mentioning that, following the 
two-fluid model 0, 0] and using the Gorter-Casimir 
temperature dependence ^4^ the limit 6(T)—> is ob- 
tained for T — > 0, as e.g. described in Ref. [3. 




FIG. 2: The lifetime r|°/r^ (n th + 1) according Eq. CT . 
This fraction only depends on the product kz and the atomic 
spin orientation. Upper figure: The spin-orientation is the 
same as in Refs.0 13, i-e. \(f\S y \i)\ 2 = |(/|5* z |i>| 2 and 
(f\S x \i) = 0. Lower figure: The atomic spin is oriented per- 
pendicular to the slab, i.e. {f\S x \i) = (f\S y \i) = 0. The 
spin-flip rate is completely suppressed for kz = in this case. 



B. Electric Dipole Transition 



The only difference between the integrals in Eqs. (|I6I17() 
and Eqs. I|8I9|) is the position of the scattering coeffi- 
cients CN(q) and Cm(^). Hence, the correction to the 
vacuum dipole-flip rate corresponding to electric field 
fluctuations for the two limits as mentioned above is 
in general opposite in sign as compared to that of the 
magnetic spin-flip case. This was also pointed out in 
Ref. 8]. It can be understood in physical terms, since the 
electric and magnetic field are perpendicular. Hence, if 
8(T) < X L (T),H,X,z or X L (T) < 8(T),H,X then the 
total electric dipole-flip rate for an atom above a slab is 
given by 
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This equation is consistent with the results in Ref . [T^ • 
Let us again consider the near-field case kz <C 1. The 
electric dipole-flip rate is then reduced to 
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provided that the atomic dipole is oriented parallel to 
the slab, i.e. provided that (/|d z |z) = 0. This result 
implies that, for an atom at the surface of the slab, i.e. 
2 = 0, there is no electric dipole-flip at all. The par- 
ticular atomic dipole orientation under consideration is 
the only one that can give zero dipole-flip rate despite the 
presence of electric field fluctuations. Furthermore, when 
the atomic spin is oriented perpendicular to the slab, the 
electric dipole-flip rate is 



r pc, 



2T° (nth + 1), 



(28) 



for the near-field case h « 1. This result was pointed 
out by Babiker in Ref. JTsj . 

To summarize, in the present paper we have reported 
results on the magnetic as well as electric decay proper- 
ties of a neutral two-level atom trapped in the vicinity of 
a dielectric body. For a slab with vacuum on both sides 
(see Fig.^l, we have obtained the flip rate for both of 
these types of transitions, for any spin or dipole orienta- 
tion. The expression for the electric and magnetic tran- 
sition rate can be solved exactly in two limiting cases, 
i.e. in the small skin depth limit for normal conducting 
metals and in the small London length limit for super- 
conductors. In these limits, the correction to the vacuum 
rate for an electric dipole transition is opposite in sign as 
compared to that of a magnetic spin transition. These 
results are consistent with well known results, e.g. Refs. 
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